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Abstract. We investigate the computational difficulty of approximating the 
partition function of the ferromagnetic Ising model on a regular matroid. Jer- 
rum and Sinclair have shown that there is a fully polynomial randomised ap- 
proximation scheme (FPRAS) for the class of graphic matroids. On the other 
hand, the authors have previously shown, subject to a complexity-theoretic as- 
sumption, that there is no FPRAS for the class of binary matroids, which is a 
proper superset of the class of graphic matroids. In order to map out the re- 
gion where approximation is feasible, we focus on the class of regular matroids, 
an important class of matroids which properly includes the class of graphic ma- 
troids, and is properly included in the class of binary matroids. Using Seymour's 
decomposition theorem, we give an FPRAS for the class of regular matroids. 

I. Introduction 

Classically, the Potts model [IT] in statistical physics is defined on a graph. Let 
q be a positive integer and G = (V, E) a graph with edge weights 7 = : e G E}; 
the weight 7 e > — 1 represents a "strength of interaction" along edge e. The g-state 
Potts partition function specified by this weighted graph is 

(1) Z Potts (G;qn)= II (l+7e*MiO,<7(i;))), 

<r:V-*[q] e={u,v}&E 

where [q] — {1, . . . , q} is a set of q spins or colours, and 5(s, s') is 1 if s = s', and 
otherwise. The partition function is a sum over "configurations" o which assign 
spins to vertices in all possible ways. We are concerned with the computational 
complexity of approximately evaluating the partition function ([Q) and generalisa- 
tions of it. For reasons that will become apparent shortly, we shall be concentrating 
on the case q — 2, which is the familiar Ising model. In this special case, the two 
spins correspond to two possible magnetisations at a vertex (or "site"), and the 
edges (or "bonds") model interactions between sites. In the ferromagnetic case, 
when 7 e > 0, for all e G E, the configurations a with many adjacent like spins 
make a greater contribution to the partition function Zp otts (G; q, 7) than those 
with few; in the antiferromagnetic case, when 7 e < 0, the opposite is the case. 
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An equivalent way of looking at (Tjrj) is as a restriction of the (multivariate) Tutte 
polynomial, which is defined as follows: 

(2) Z(G;q, 1 ) = J27Aq K(A) - lVl , 

ACE 

where = HeeA^e anc ^ denotes the number of connected components in 

the graph (V, A) [T31 (1-2)]- This is perhaps not the most usual expression for the 
multivariate Tutte polynomial of a graph, but it conveniently generalises to the 
Tutte polynomial of a matroid fl5]), also [131 (1-3)], which is the main subject of 
this article. Although (J2J) and (0Q) are formally quite different, they agree when q 
is a positive integer, up to a factor of q~\ y \. 

Jaeger, Vertigan and Welsh [7] were the first to consider the computational 
complexity of computing the Tutte polynomial. They considered the classical 
bivariate Tutte polynomial in which the edge weights are constant, i.e., 7 e = 7 
for all e G E. Their approach was to fix q and 7, and consider the computational 
complexity of computing ([2]) as a function of the instance graph G. Jaeger et al. 
showed, amongst other things, that computing the Tutte polynomial exactly is 
#P-hard when q > 1 and 7 G (—1, 00) — {0}. In particular, this means that the 
partition function (CQ) of the Potts model is computationally intractable, unless 
P = #P. 

In the light of this intractability result, it is natural to consider the complexity of 
approximate computation in the sense of "fully polynomial approximation scheme" 
or FPRAS. Before stating the known results, we quickly define the relevant con- 
cepts. A randomised approximation scheme is an algorithm for approximately 
computing the value of a function / : E* -> 1. The approximation scheme has a 
parameter e > which specifies the error tolerance. A randomised approximation 
scheme for / is a randomised algorithm that takes as input an instance x G X* 
(e.g., for the problem of computing the bivariate Tutte polynomial of a graph, the 
graph G) and a rational error tolerance e > 0, and outputs a rational number z (a 
random variable of the "coin tosses" made by the algorithm) such that, for every 
instance x, 

(3) Pr [e~ £ f(x) < z < e £ f(x)} > 1 . 

The randomised approximation scheme is said to be a fully polynomial randomised 
approximation scheme, or FPRAS, if it runs in time bounded by a polynomial in 
\x\ and Note that the quantity 3/4 in Equation ()3|) could be changed to any 
value in the open interval (|, 1) without changing the set of problems that have 
randomised approximation schemes [9], Lemma 6.1]. 

For the problem of computing the bivariate Tutte polynomial in the antiferro- 
magnetic situation, i.e., — 1 < 7 < 0, there is no FPRAS for (J2]) unless RP = 
NP [1]. This perhaps does not come as a great surprise, since, in the special case 
7 = — 1, the equivalent expression ([T]) counts proper colourings of a graph. So 
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we are led to consider the ferromagnetic case, 7 > 0. Even here, Goldberg and 
Jerrum [5] have recently provided evidence of computational intractability (under 
a complexity theoretic assumption that is stronger than RP 7^ NP) when q > 2. 

The sequence of results so far described suggest we should focus on the special 
case q = 2 and 7 > 0, i.e., the ferromagnetic Ising model. Here, at last, there is a 
positive result to report, as Jerrum and Sinclair [8] have presented an FPRAS for 
the partition function ([2]), with q = 2 and arbitrary positive weights 7. As hinted 
earlier, the Tutte polynomial makes perfect sense in the much wider context of 
an arbitrary matroid (see Sections [2] and [3] for a quick survey of matroid basics, 
and of the Tutte polynomial in the context of matroids). The Tutte polynomial of 
a graph is merely the special case where the matroid is restricted to be graphic. 
It is natural to ask whether the positive result of [3J extends to a wider class of 
matroids than graphic. One extension of graphic matroids is to the class of binary 
matroids. Goldberg and Jerrum [5] recently provided evidence of computational 
intractability of the ferromagnetic Ising model on binary matroids, under the same 
strong complexity-theoretic assumption mentioned earlier. 

Sandwiched between the graphic (computationally easy) and binary matroids 
(apparently computationally hard) is the class of regular matroids. Since it is 
interesting to locate the exact boundary of tractability, we consider here the com- 
putational complexity of estimating (in the FPRAS sense) the partition function 
of the Ising model on a regular matroid. We show that there is an FRPAS in this 
situation. (See Section [2] for matroid definitions and Section [3] for the definition of 
the Tutte polynomial Z(Ai,-y) of a matroid Ai.) 

Theorem 1. There is an FPRAS for the following problem. 

Instance: A binary matrix representing a regular matroid At. A set 7 = 

{7 e : e G E(Ai)} of non-negative rational edge weights. 
Output: Z(M;2,j). 

Aside from the existing FPRAS for graphic matroids, which also works, by du- 
ality, for so-called cographic matroids, the main ingredient in our algorithm is 
Seymour's decomposition theorem for regular matroids. This theorem has been 
applied on at least one previous occasion to the design of a polynomial-time algo- 
rithm. Golynski and Horton [5] use the approach in their algorithm for finding a 
minimum-weight basis of the cycle space (or circuit space) of a regular matroid. 
The decomposition theorem states that every regular matroid is either graphic, 
cographic, a special matroid on 10 elements named Rio, or can be decomposed as 
a certain kind of sum (called "1-sum", "2-sum" or "3-sum") of two smaller regular 
matroids (see Theorem [3]). Since we know how to handle the base cases (graphic, 
cographic and Rio), it seems likely that the decomposition theorem will yield a 
polynomial time algorithm quite directly. However, there is a catch (which does 
not arise in j6]). When we pull apart a regular matroid into two smaller ones, 
say into the 3-sum of Ati and AA2, four subproblems are generated for each of the 
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parts Ai\ and A^- This is fine if the decomposition is fairly balanced at each step, 
but that is not always the case. In the case of a highly unbalanced decomposition, 
we face combinatorial explosion. 

The solution we adopt is to "solve" recursively only the smaller subproblem, 
say A^2- Then we construct a constant size matroid I3 that is equivalent to AA2 in 
the context of any 3-sum. We then glue X 3 onto Aii, using the 3-sum operation, 
in place of Ai2- The matroid X3 is small, just six elements, and has the property 
that forming the 3-sum with Ai\ leaves Ai\ unchanged as a matroid, though it 
acquires some new weights from X 3 . (In a sense, X 3 is an identity for the 3-sum 
operation.) Then we just have to find the partition function of Ai\ (with amended 
weights). Since we have four recursive calls on "small" subproblems, but only one 
on a "large" one, we achieve polynomially bounded running time. 

The fact that X3 is able to simulate the behaviour of an arbitrary regular matroid 
in the context of a 3-sum is a fortunate accident of the specialisation to q — 2. 
Since the existing algorithm for the graphic case also makes very particular use of 
the fact that q = 2, one gets the impression that the Ising model has very special 
properties compared with the Potts model in general. 

2. Matroid preliminaries 

A matroid is a combinatorial structure that has a number of equivalent defi- 
nitions, but the one in terms of a rank function is the most natural here. A set 
E (the "ground set") together with a rank function r : E — > N is said to be 
a matroid if the following conditions are satisfied for all subsets A, B C E: 
(i) < r(A) < \A\, (ii) A C B implies r(A) < r(B) (monotonicity) , and 
(iii) r(A U B) + r(A D B) < r(A) + r(B) (submodularity). A subset ACE 
satisfying r(A) = \A\ is said to be independent; a maximal (with respect to in- 
clusion) independent set is a basis, and a minimal dependent set is a circuit. 
A circuit with one element is a loop. We denote the ground set of matroid Ai 
by E(Ai) and its rank function by tm- To every matroid Ai there is a dual ma- 
troid Ai* with the same ground set E = E(Ai) but rank function rjw* given by 
rM*(A) = \A\ + tm{E — A) — tm(E). A cocircuit in Ai is a set that is a circuit 
in equivalently, a cocircuit is a minimal set that intersects every basis. A co- 
circuit with one element is a coloop. A thorough exposition of the fundamentals 
(and beyond) of matroid theory can be found in Oxley's book [TU] . 

We give names to two specific matroids which we will use later. A/i is defined to 
be the matroid with ground set E{Af\) = {p} in which r^ip) = 1. A3 is defined to 
be the matroid with ground set E(Afs) = {pi,P2,P3} in which £7 (A/3) is a circuit. 

Important operations on matroids include contraction and deletion. Suppose 
T C E is any subset of the ground set of matroid Ai. The contraction Ai/T 
of T from Ai is the matroid on ground set E — T with rank function given by 
r M/T (A) = r M (A[jT)-r M (T), for all A C E-T. The deletion M\T of T from M 
is the matroid on ground set E — T with rank function given by rj^/j-(A) = rj^(A), 
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for all A C E — T. These operations are often combined, and we write M/T\S 
for the matroid obtained by contracting T from M and then deleting S from 
the result. The operations of contraction and deletion are dual in the sense that 
(M\T)* = M* /T '. For compactness, we shall often miss out set brackets, writing 
M/pi\p 2 ,P3, for example, in place of M/{pi}\{p 2 ,ps}. The restriction M \ S 
of M to S C E is the matroid on ground set S that inherits its rank function 
from M; another way of expressing this is to say M \ S — M\{E — S). 

The matroid axioms are intended to abstract the notion of linear independence of 
vectors. Some matroids can be represented concretely as a matrix M with entries 
from a field K, the columns of the matrix being identified with the elements of E. 
The rank r(A) of a subset A C E is then just the rank of the submatrix of M 
formed from columns picked out by A. A matroid that can be specified in this way 
is said to be representable over K. A matroid that is representable over GF(2) 
is binary, and one that is representable over every field is regular; the regular 
matroids form a proper subclass of binary matroids. Another important class of 
matroids are ones that arise as the cycle matroid of an undirected (multi)graph 
G = (V, E). Here, the edge set E of the graph forms the ground set of the matroid, 
and the rank of a subset A C E is defined to be r(A) — \V\ — k(A), where k(A) is 
the number of connected components in the subgraph (V, A). A matroid is graphic 
if it arises as the cycle matroid of some graph, and cographic if its dual is graphic. 
Both the class of graphic matroids and the class of cographic matroids are strictly 
contained in the class of regular matroids. 

Now for some definitions more specific to the work in this article. A cycle in 
a matroid is any subset of ground set that can be expressed as a disjoint union 
of circuits. We let C(M) denote the set of cycles of a matroid M. If M is a 
binary matroid, the symmetric difference of any two cycles is again a cycle [TO] 
Thm. 9.1.2(vi)], so C(M), viewed as a set of characteristic vectors on E(M), 
forms a vector space over GF(2), which we refer to as the circuit space of M. 
Indeed, any vector space generated by a set of vectors in GF(2) E can be regarded 
as the circuit space of a binary matroid on ground set E (see the remark following 
[TP] Cor. 9.2.3]). The set of cycles C(M) of a matroid M determines the set 
of circuits of M (these being just the minimal non-empty cycles), which in turn 
determines M. Note, from the definition of "cycle", that 

C{M\T) = {C e C(M) | C c E\T}. 

The term "cycle" in this context is not widespread, but is used by Seymour [12] 
in his work on matroid decomposition; the term "circuit space" is more standard. 

Consider two binary matroids Mi and M 2 with E(Aii) — EiUT and E(Ai 2 ) = 
E 2 UT and E x D E 2 = 0. The delta-sum of Mi and M 2 is the matroid Mi A M 2 
on ground set Ei U E 2 with the following circuit space: 

C(MiAM 2 ) = {C C EiUE 2 : C = C X ®C 2 for some G x G C(Mi) and C 2 G C{M 2 )}, 
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where © denotes symmetric difference [12]. (The right-hand side of the above 
equation defines a vector space over GF(2), and hence does describe the circuit 
space of some matroid.) A situation of particular interest occurs when Mi \ T = 
M2 | T = N (say), i.e., Mi and M 2 have the same restriction, namely N, to T. 
To investigate this situation, consider the following collection of cycles, this time 
on the larger set E\ U T U E 2 : 
(4) 

C(A) = {C C E 1 UTUE 2 : C = Ci®C 2 for some d G C(Mi) and C 2 G C(M 2 )}, 

For the same reason as before, this is the circuit space of a matroid, which we'll 
denote by A. The following lemma shows that A is an amalgam of .Mi and M2] 
that is to say A | (E 1 U T) = Mi and A | (i? 2 U T) = M 2 . Note that there is no 
reason a priori to suppose that an amalgam exists, though the lemma shows that 
it does in this case. 

Lemma 2. Suppose Mi and M 2 are binary matroids with Mi \ T = M 2 \ T, 
where T = E(Mi) PI E(M 2 ). The the matroid A with circuit space given by 
(with Ei = E(Mi) — T and E 2 = E(A4 2 ) — T) is an amalgam of M.\ and M. 2 . 

Proof. Certainly, any cycle in M.\ is also a cycle in A \ (E\ U T). Now consider 
any cycle in A \ (E\ U T); this cycle has the form C\ © C 2 where C\ G C(Aii), 
C 2 G C(M. 2 ) and C 2 C T. But since M. l \ T = M. 2 \ T, we see that C 2 is also a 
cycle in Mi, and hence C = C 1 © C 2 is a cycle in M x . Thus, C(A \ (E x U T)) = 
C(M\), and A | (Ex U T) — M\. By symmetry, a similar identity holds for M 2 , 
and we deduce that A is an amalgam of Mi and M 2 . □ 

There is a routine method that gives an upper bound on the rank function of 
the amalgam A in terms of the constituent matroids M\ and M 2 . Since 

C(A\T) = {C e C{A) | C C £i U E 2 } = C(Mi A M 2 ), 

we know Mi A .M2 = A\T. So the rank function of A immediately gives us a 
handle on the rank function of the delta-sum Mi A M 2 . We shall make use of 
this computational tool in the following section. 

The special case of the delta-sum when T = is called the 1-sum of Mi and M 2 - 
(It is just the direct sum of the two matroids.) The special case when T is a 
singleton that is not a loop or coloop in Mi or M 2 , and |J5(.Mi)|, \E(M 2 )\ > 3, 
is called the 2-sum of Mi or M 2 and is denoted Mi ©2 M 2 . Finally, the special 
case when |T| = 3, T is a circuit in both Mi and M 2 but contains no co-circuit 
of either, and \E(Mi)\, \E(M2)\ > 7, is called the 3-sum of Mi or M 2 and is 
denoted Mi@2,M 2 p2]. 

Our main tool is the following celebrated result of Seymour [121 (14.3)]: 

Theorem 3. Every regular matroid M may be constructed by means of 1- 2- and 
3-sums, starting with matroids each isomorphic to a minor of M, and each either 
graphic, cographic or isomorphic to a certain 10-element matroid Riq. 
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It is important for us that a polynomial-time algorithm exists for finding the de- 
composition promised by this theorem. As Truemper notes [13] , such an algorithm 
is given implicitly in Seymour's paper. However, much more efficient algorithms 
are known. In particular, [T4"| (10.6.1)] gives a polynomial-time algorithm for test- 
ing whether a binary matroid is graphic or cographic. Also, [15] gives a cubic 
algorithm for expressing any regular matroid that is not graphic, cographic or 
isomorphic to Rio as a 1-sum, 2-sum or 3-sum. For an exposition of this result, 
[T4J (8.4.1)] gives a polynomial-time algorithm that takes as input a matrix repre- 
senting a binary matroid and produces a 1-sum or a 2-sum decomposition, or, if 
the matroid is 3-connnected, produces a 3-sum decomposition if one exists. From 
the proof of [321 (14.3)], and the conditions on the matroid, (that it be regular, 
but not graphic, cographic or isomorphic to -Rio), the matroid has a 3-separation 
in which the parts are sufficiently large, so by [HI (8.3.12)], the desired 3-sum 
decomposition does exist (so is constructed). 

Truemper's definitions are slightly different from Seymour's (which we use). 
However, it is easy to check that Truemper's 1-sum [T4"| Section 8.2] is a 1-sum in 
the sense of Seymour. It can also be checked that Truemper's 2-sum is a 2-sum 
in the sense of Seymour. (For this it helps to note that the ground-set element in 
E(Aii) n E(Ai 2 ) is represented by the column corresponding to element u x" in 
Truemper's matrix / | B 1 and by the column corresponding to element u y" in his 
matrix I | B 2 .) Truemper's 3-sums [T4"| Section 8.2] are not quite 3-sums in the 
sense of Seymour, but, as Golynski and Horton have noted [6], it is easy to apply 
an exchange operation [T3J Section 8.5] (in linear time) to obtain a 3-sum in the 
sense of Seymour. 

While Seymour's decomposition theorem is in terms of 1-sums, 2-sums and 3- 
sums, it will be convenient for us to do our preparatory work, in the following 
section, using the slightly more general notion of a delta-sum. 

3. TUTTE POLYNOMIAL AND DECOMPOSITION 

Suppose M. is a matroid, q is an indeterminate, and 7 = : e G E(M)} is 
a collection of indeterminates, indexed by elements of the ground set of M.. The 
(multivariate) Tutte polynomial of M. and 7) is defined to be 

(5) Z(M;q,~f)= X ™~™ {A) > 

ACE(M) 

where 7a = FLeA^e |13> (1-3)]- in this article, we are interested in the Ising model, 
which corresponds to the specialisation of the above polynomial to q = 2, so we 
shall usually omit the parameter q in the above notation and assume that q is set 
to 2. As we are concerned with approximate computation, we shall invariably be 
working in a environment each of the indeterminates 7 e is assigned some real value; 
furthermore, we shall be focusing on the ferromagnetic case, in which those values 
or weights are all non- negative. For convenience, we refer to the pair (Ai,-y) as a 
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"weighted matroid", and we will assume throughout that 7 e > for all e G E(M). 
For background material on the multivariate Tutte polynomial, and its relation to 
the classical 2- variable Tutte polynomial, refer to Sokal's expository article [13] . 

In order to exploit Theorem [3] we need to investigate how definition (J5]) behaves 
when M = Mi A M 2 . 

Suppose that Ex, T and E 2 are mutually disjoint sets and consider two binary 
matroids Mx and M 2 with E(Mx) = E x U T and E(M 2 ) = E 2 UT. Assume, 
as we did in Section [2] that M.x and Ai 2 have the same common restriction 
iV = A4i | T = Ai 2 | T to T. Then, as we saw in Section [2] the delta-sum 
M. — M.i A M. 2 may be expressed as M. = A\T where A is an amalgam of Aix 
and M 2 . 

Oxley [ini Section 12.4] has given an upper bound on the rank function of any 
amalgam of .Mi and Ai 2 . Thus, the upper bound applies to r^, and hence tm- 
In particular, 

(6) VIC^U E 2 ,r M (X) = r A (X) < mm{r](Y) :YDX}, 
where 

V (Y) = r Ml (Y n S(Mi)) + r M2 (Y n S(X 2 )) - r^(F n T). 
For i C £ ; and S 1 C T, let e 4 (A, S) = r Mi (A US)- r Mi (A). Intuitively, this 
is the "excess" rank that S adds to A. For A C Ex and B C _E 2 , let c(A, S) = 
r M (AUB) — r Ml (A) — r M , 2 (B). Intuitively, c(A, B) is the appropriate "correction" 
to the rank function of M. relative to the rank function of the simple 1-sum of Aii 
and M 2 . 

Equation ()6]) gives 

r M {A UB)< min (r Ml (A U S) + r M2 (B US)- r N (S)) 

= r Ml (A) +r M2 (B) + min (ex(A, S) + e 2 (B, S) - r N (S)) , 
so Equation (E]) gives 

(7) c{A, B) < min {ex(A, S) + e 2 (B, S) - r N (S)) . 

From the definitions of the matroids Mx and A3, we can now see that, when 
N = ATx, 

(8) c{A,B) < min(0,ex(A,p) + e 2 (B,p) - 1). 
Also, when N = A/3, 

(9) c(A, B) < min (ex(A, S) + e 2 (B, S) - \S\) . 

S dT 

Here we use the fact that the right-hand side of ([7]) cannot have a unique minimum 
at S = T, in order to restrict the minimisation to strict subsets of T.) 

The interaction between the classical bivariate Tutte polynomial and 2- and 3- 
sums has been investigated previously, for example by Andrzejak [I]. In principle, 
it would be possible to assure oneself that his proof carries over from the bivariate 
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to the multivariate situation, and then recover the appropriate formulas (Lemmas 
H] and below) by appropriate algebraic translations. However, in the case of the 
3-sum there is an obstacle. While Andrzejak's identities remain valid, as identities 
of rational functions on q, they become degenerate (through division by zero) under 
the specialisation q = 2. In fact, this degeneracy is crucial to us, in that it reduces 
the dimension of the bilinear form in Lemma |5]from five, as in Andrzejak's general 
result, to four, in our special formula for q = 2. In light of these considerations, and 
because certain intermediate results in our proofs are in any case required later, 
we derive the required formulas here. First the formula for (a slight relaxation of) 
the 2-sum. 

Lemma 4. Suppose (A^i,^) and (A^2>72) are weighted binary matroids with 
E(J\4i) Pi E(M. 2 ) = {p}, where p is not a loop in either M.\ or M. 2 . Denote by 
7 = 7j A 7 2 the weighting on E{M.\ AM 2 ) inherited from A^i and M. 2 . Then 
(10) 

Proof. The machinery of the proof is a little heavy in relation to the scale of the 
result, but its use will provide a warm-up for the proof of the analogous result for 
3-sums (Lemma [5]). 

Write E{M\) as E\ U {p} and E(M 2 ) as E 2 U {p} and note that 

N = M 1 \{p}=M 2 \ {p}=M- 
Use the definition of c(A, B) to write the Tutte polynomial of the delta-sum as 

(11) Z{Mi AM 2 ;j) = Yj^AlBq~ rMl{A) ~ rM * {B) q~ <A ' B) - 

ACE! BcE 2 

We have from flH]), that 

'-1 iff e 1 (A,{p}) = and r 2 ( B.{ P }) = 0: 
otherwise. 



c{A,B) < 



We show in Appendix 16.1.11 that this upper bound is also a lower bound, so 
c(A,B) 



-1 iSe 1 (A,{p}) = 0ande 2 (B,{p}) = 0; 
otherwise. 



Thus, letting z be the column vector z l = (zq, z\) T , for i e {1,2}, where 
zh= J2 lAq~ rM > {A) and z\= ^ l A q~ r ^ {A \ 

A:ACEi,ei(Mp})=0 A:ACEi, ei {A{p})=^ 

we may express identity (11 II) as 

(12) Z(M 1 AM 2 ;-f) = (z 1 ) T Cz 2 , 
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Having expressed the left-hand side of (flOj) in terms of z 1 and z 2 , we now wish 
to do the same for the right-hand side. For % e {1, 2}, let 

z* = (4,4) T = [Z {Mi\p; 7J, Z(Mi/p; 7i)) T , 

and observe that 



ACE,; 



and 



so 



(13) 
where 



(We are using the fact that 7 agrees with ~f 1 on and with 7 2 on £7 2 -) Now, 
Zq > and 5} > (since the weights 7 e are non-negative and the term A = 
contributes positively to the sum defining z\). Together with ([TBI these facts 
imply 

(14) 4<z[< 24 

an inequality we shall need later. With 

D = 

we can express the right-hand side of (flQj) . using transformation (LTBl . as 



(15) 



V 1\T 



£> z 2 = (F^Wi 



^l^T 



(V T W)5 2 . 



Comparing (|T2|) and ffTBI) . we see that to complete the proof, we just need to verify 
the matrix equation C = V T DV with C and V specialised to q — 2, and this is 
easily done. 

The same calculation can be carried through for 
3-sum. Denote by D the matrix 

/ 4 -1-1 -l\ 
-110 



□ 

(a slight relaxation of) the 



D 



V 








1 
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Lemma 5. Suppose that (.Mi, 7l ) and (M2, I2) are weighted binary matroids with 
E(Aii) PI E(M2) = T = {pi,P2,Ps} , and suppose also that T is a circuit in both 
Mi and M 2 - For i = 1,2, let 

z l = {Z(M l \T-j i ),Z(M i /p 1 \p2,P3]li),Z(Mi/p2\pi,p 3 ]j i ),Z(M i /p 3 \p 1 ,p 2 -,j i )) T . 

Denote by 7 = / y 1 A 7 2 the weighting on E(M\ A M2) inherited from Mi and 
M 2 - Then 

(16) Z(Mi A M 2 ] 7) = (z x ) T Dz 2 . 

Proof. Note that N = Mi \ T = M 2 \ T = J\f 3 . Use the definition of c(A,B) to 
write the Tutte polynomial of the delta-sum as 

(17) Z(M 1 AM 2 ;~f)=J2 iAlBq~ rM ^ A) ~™* {B) q~ c{A,B) ■ 

AC_Bi BcE 2 

Now, to construct an equation analogous to ffTUl) . we will work with the minors 
Mi/Si\S 2 for partitions (Si, S2) of T. This gives eight different minors of Mi, but 
some of them are equivalent. In particular, since T is a circuit of Mi (which means 
that pi, P2 and ^3 are linearly dependent in the binary matrix representing Mi 
but any pair of them is independent), any partition (Si, S2) with \Si\ > 2 gives a 
matroid Mi/Si\S2 which is equivalent to Mi/T. To see this, note that for X C Ei, 
r Mi / Sl \s 2 (X) = r Mi (XUS 1 )-r Mi (S 1 ) = r Mi (X U T) - 2 = r Mi (X U T) - r Mi (T) . 
Thus, we will use the five minors with |Si| 7^ 2. 

Let's collect the formulas for the Tutte polynomials of the minors that we'll 
need. 

(18) Z(M\T; 7 *)= E^" r " l(A) 

ACE, 

Z(Mi/ Pj \T - Pf, id = lA q' rM ^ A) q l - ei{MPi} \ for 3 G {1, 2, 3}, 

ACE, 

and 

Z(M i /T ] y i )= J2lAq- r ^ iA) q 2 - edA ' T) 

ACE, 

With an eye to the proof of Lemma IU we next need to understand the function 
ei(A,S). It turns out that, with A fixed, ei(A,S) is completely determined by 
its value on the on the singleton sets S = {pi}, {P2}, {^3}- Not only that, there 
are just five possible for the values of ei(A,S) on those singletons. For i = 1,2, 
the following collection of predicates {ip\, . . . , (p\} on A C E i} captures those five 
possibilities. 

PO. cpl(A) iff ei(A, {pi}) = d(A, {p 2 }) = e t (A, {p 3 }) = (in which case, e^A, S) = 
for all S). 

PI. <p\(A) iff ei(A,{pi}) = and e^A, {p 2 }) = e { (A, {p 3 }) = 1 (in which case, 
ei(A,S) = 1 when \S\ > 2). 
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P2. (fi(A) iff ei(A,{p 2 }) = and e^A, { Pl }) = ei (A,{p 3 }) = 1 (in which case, 

ei(A,S) = 1 when \S\ > 2). 
P3. <p\(A) iff ej(v4, {p 3 }) = and ei(A, {pi}) = ei(A,{p 2 }) = 1 (in which case, 

ei(A,S) = 1 when \S\ > 2). 
P4. ^(A) iff ei{A, {pi}) = ei(A, {p 2 }) = e;(A, {p 3 }) = 1 (in which case, e^A, S) = 

2 when |S| > 2). 

Of course, in all cases, Ci(A, 0) = 0. 

First observe that it is not possible that exactly one of e^A, {pi}), e^A, {P2}) 
and ei(A, {p 3 }) to take the value 1, since this would imply that two members of T 
are dependent on A, and hence the third would be. Thus, FtD]-F|l] are exhaustive 
as well as mutually exclusive. 

We still need to verify the additional information provided in parentheses. In 
case PEJ &i{A, S) = for all S, since all of the elements in T are dependent 
on A. In case FTf] (F(2] and Fj3] are symmetrical) p± is is dependent on A, and 
so e i (A,{p 1 ,p 2 }) = 6i(A, {p 2 }) = 1; then e<(A,T) = e i (A,{p 1 ,p 2 }) = 1, since p 3 
depends on pi and p 2 . The final parenthetical claim (in PH]) is slightly trickier, and 
relies on the fact that the matroids we are working with are binary. (There is a 
simple counterexample based on the uniform matroid U^ 2 for the claim in general.) 
So consider the elements of A as columns of the representing matrix (vectors over 
GF(2)), and similarly consider the elements of S as columns. We must rule out 
the possibility that Ci(A, S) = 1 for some 2-element subset S, say S = {pi,P2}- So 
suppose ei(A,{pi}) = ei(A,{p 2 }) = e i (A,{p 1 ,p 2 }) = 1. Then p 2 is dependent on 
A U {pi}, i.e., so (also viewing p\ and p 2 as vectors over GF(2)), p 2 = p\ + YleeB e ' 
for some B C A. (Note that p\ must be included in this expression, otherwise p 2 
would be dependent on A.) Since T is a circuit we know that p 3 = p\ + p 2 . Now 
write p 3 = pi + p 2 = pi + pi + J2eeB e = Sees e - Thus p 3 is dependent on A, and 
ei(A,p 3 ) = 0, contrary to what we assumed. 

By analogy with the proof of Lemma H] define z % = (z l Q , . . . , z\) T for i = 1, 2 by 

4 = E 7A?- r ^ (A) 

and 

z l = (Z(M\T;7i),^(MM\P2,P3;7i), 

Z(M l /p 2 \ P i,p 3 ; 7j ), Z(Mi/p 3 \p u P2; 7,), ^(^,/T; 7 ,)) T - 

Note that the dimension here is one greater (five instead of four) than in the 
statement of the lemma, but this extra dimension will disappear towards the end 
of the proof. As before, we can re-express the left-hand side of ([IB]) in terms of 

z\i\ 

Our first step will be to show that the value c(A, B) depends on the predicates 
<Pj(A) and tp\{B) but does not depend otherwise on A and B. Thus, we will be 
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able to define a 5 x 5 matrix C with Cj^ = q~ c ( A < B \ for < j, k < 4, where A C E\ 
and B C E 2 are any sets satisfying y?}^) and ip\(B). Using Equation (j^J), we can 
construct an element-wise lower bound for the matrix C . In particular, 



(q 2 



c > 



q q q 



q q 

q 
l 



l l 

1 q 1 

1 1 q 

1 1 1 



1\ 
1 
1 
1 



For example, if ip\(A) and <p 2 {B) are true, then ei(A, {pi}) = and ei(A, {p 2 }) = 
e 1 {A,{p 3 }) = 1 so ei(A,S) = 1 when \S\ > 2. Also, e 2 (S, {p 2 }) = and 
e 2 (S,{pi}) = e 2 (B,{p 3 }) = 1 so e 2 {B,S) = 1 when |5| > 2. Now Equa- 
tion says c(A,i?) < min^c-r (e\(A, S) + e 2 (B, S) — \S\). Checking the pos- 
sibilities for S, we find c(A, B) < 0, so g~ c ( A ' B ) > 1. Thus, the entry in row 1 (the 
second row from the top) and column 2 (the third column from the left) above 
is 1. The rest of the entries are filled in in the same way. 

We show in Appendix 16.1.21 that this lower bound is also an upper bound, so 



(q 2 



c 



q q q 



q q 
q 
q 
l 



l l 

1 q 1 

1 1 q 

1 1 1 



1\ 

1 
1 
1 



(A 2 2 2 l\ 

2 2 111 

2 12 11 

2 112 1 

\1 1 1 1 1/ 



Then the left-hand side of (|T6l) may be expressed as 



(19) Z(M 1 AM 2 ;j) = (z 1 ) L Cz 2 . 

Also as before, there is a linear relationship between 



and z\ expressed by 



(20) 
where 



(1 


l 


l 


l 


i\ 




A 


1 


1 


1 


A 


q 


q 


l 


l 


l 




2 


2 


1 


1 


i 


q 


l 


q 


l 


l 




2 


1 


2 


1 


i 


q 


l 


l 


q 


l 




2 


1 


1 


2 


i 


\q 2 


q 


q 


q 


V 




V 4 


2 


2 


2 


V 



V 



It is straightforward to verify correctness of V row by row. For example, the third 
row expresses the easily checked fact that 

Z(M i /p 2 \p 1 ,p 3 ; 7j = qzl + z[ + qz\ + z\ + 5j. 

Identity fl20|) allows us to rewrite the right-hand side of f|T6l) in terms of z 1 and z 2 . 
Let D' be the 5x5 matrix obtained from D in the statement of the theorem by 
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padding it out to the right and below by a single column and row of zeros. Then 
the right-hand side of f[T5]) may be written as 

(21) (z'fDz 2 = (Vz'fD'Vz 2 = (z'f^D'V)^. 

Comparing ( TI9]) and fl2T|) . we see that to complete the proof, we just need to verify 
the matrix equation 

(22) C = V T D'V, 

with C and V specialised to q — 2, and this is easily done. □ 

Remark 1. With q regarded as an indeterminate, Equation 1(22]) has a unique 
solution in which D' is a matrix of full rank, whose entries are rational functions 
of q. When q is specialised to 2, the rank of matrix C drops to 4, allowing a 
solution in which D' also has rank 4- When q = 2, there is flexibility in the choice 
of the matrix D' , which we exploit in order to reduce the dimension by 1. This is 
the reason that the right-hand side of (TT^j is of dimension 4, whereas one would 
expect dimension 5 a priori. This apparent accident is crucial to the design of the 
algorithm. 

4. Signatures 

The goal of this section is to show that, for every weighted binary weighted 
matroid (Ai, 7) with distinguished element p, there is a small (2-element) weighted 
matroid (X 2 , 8) that is equivalent to M. in the following sense: if we replace (.M, 7) 
by (X 2 , S) in the context of any 2-sum, the Tutte polynomial (specialised to q — 2) 
of the 2-sum is changed by a factor that is independent of the context. Moreover, 
the weighted matroid (X 2 ,<5) can readily be computed given a "signature" of M.. 
There is also a 6-element weighted matroid (I3, 8) that does a similar job for 
3-sums. 

To make this precise, let X 2 be the matroid with a 2-element ground set {p, e} 
that form a 2-circuit. X 2 can be viewed as the cycle matroid of a graph consisting 
of two parallel edges, p and e. 

Lemma 6. Suppose (Ai,~f) is a weighted binary matroid with a distinguished 
element p that is not a loop. Then there is a weight d > such that, for every 
weight function 6 with 5 e = d, 

Zjl 2 /p-8) _ ZjM/p-n) 
Z(1 2 \ P] S) Z{M\p;i)' 

The value d can be computed from Z(Ai\p;~f) and Z(Ai/p;~f) - it does not 
otherwise depend upon M. or 7. 
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Proof. Let z = Z(Ai\p; 7) and Z\ = Z(A4/p; 7), and set d = 2{zi — z )/(2z — z\). 
Inequality (1141) implies that the numerator is non-negative, and the denominator 
is strictly positive. Then 

Z(T 2 /p- d)/Z(Z 2 \p; S) = (l + d)/(l + \d) = Zl /z = Z(M/p; y)/Z(M\p; 7)- 

□ 

Corollary 7. Suppose (Mi,"fi) and (.M 2 ,7 2 ) are weighted binary matroids with 
E(Aii) fl E(AA 2 ) = {p}, where p is not a loop in either Ai\ or Ai 2 . Then there 
is a weight d > such that, for every weight function 8 with 5 e = d, 

Z{M 1 A M 2 ; 7l A 7a ) = (Z{M X A Z 2 ; 7l A S), 

where ( = 2(2+d)~ 1 Z(A4 2 \p; j 2 ). The value d can be computed from Z(A4 2 \p; 7 2 ) 
and Z(Ai 2 /p;j 2 ) — it does not otherwise depend upon (.Mi, 7^ or (.M 2 ,7 2 ). 



Z(MiAM 2 ; 7iA7 2 ) = {Z{M 1 \ P]ll ),Z{M 1 /p; ll )) ( 2 , (]) Z{M 2 \p; l2 ). 



Proof. Let s = Z(M 2 /p;j 2 )/Z(M 2 \p;j 2 ). By LemmaU 

' 2 

-1 1 J \sj 
Similarly, using Lemma [6J 

Z{Mi A X 2 ; 7l A S) = (Z(Mi\p; 7l ), Z{Mi/p; 7i)) Q Q ^(X 2 \p; 5). 

Thus, C = ZCM 2 \p; 7 2 )/^(^\p; 5). □ 

It is crucial in Corollary [7] that £ does not depend at all on (-Mi, 7l ). As we 
shall see presently, an analogous result holds for 3-sums, though the calculations 
are more involved. For a weighted binary matroid (.M,7) with distinguished 
elements T = {pi,p 2 ,P3j-, the signature of (M.,~y) with respect to T is the vector 

Z(M/p 1 \p 2 ,p 3 ;j),Z(M/p 2 \pi,p 3 ;7),Z(M/p 3 \p 1 ,p 2 ;j) 

a(M;l,j) = ~ . 

Z(M\T; 7 ) 

What we are seeking in the 3-sum case is a small matroid whose signature is equal 
to that of a given binary matroid Ai. Such a matroid will be equivalent to Ai 
in the context of any 3-sum. Before constructing such a matroid we give two 
technical lemmas that investigate inequalities between the Tutte polynomial of 
various minors of a binary matroid. These inequalities will restrict the domain of 
possible signatures that can occur. 
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Lemma 8. Suppose that (Ai,~f) is a weighted binary matroid with distinguished 
elements T = {pi,p 2 ,p 3 } , and that T is a circuit in M.. Let 

Z T = (zq, Zi, Z 2 , Z 3 , Z4) 

= (Z(M\T ] -f),Z(M/ Pl \p 2 ,p 3 ;j), 

Z(M/p 2 \p llP3 ; 7), Z(M/p 3 \ Pl ,p 2 ; 7), Z(M/T; 7)) . 

Then the following (in) equalities hold: (i) z > 0, (ii) z < zi,z 2 ,z 3 < 2zo, 
(Hi) ~z± < zi, z 2 , z 3 < Zi and (iv) z 1 + z 2 + z 3 = 2z + z 4 . 

Proof. Inequality (i) follows from the definition §5§ since the contribution of A = 
is 1 and the contribution of every other A is non- negative. Consider identity (120]) . 
Since z > (coordinatewise), the vector z in the statement of the lemma is in the 
cone generated by the columns of V (i.e., z is a non- negative linear combination 
of columns). Denote the rows of V by Vo, . . . , V4. Then, interpreting all vector 
inequalities coordinatewise, (ii) is a consequence of Vo < V\, V 2 , V 3 < 2Vo (iii), with 
non-strict inequality, of ^V^ < VI, V 2 , V 3 < V4, and (iv) of Vi + V 2 + V 3 = 2V + V4. 
We know that z^ > 0, since the weights are non- negative and its expansion contains 
at least one non- zero term, namely the one corresponding to A = 0; this implies 
that the first inequality in (iii) must be strict. □ 

Lemma 9. Under the same assumptions as Lemma® ?iz 2 , z±z 3 , z 2 z 3 < ZqZ±. 

Proof. Let the ground set of A4 be E U T where where E fl T = 0. Recall the 
earlier definition 

z k = lA<T rM{A \ for < k < 4, 

A:AeE,<p k (A) 

and recall also that the the predicates <fk are exhaustive and mutually exclusive, 

so that Z(M\T; 7) = z + z\ H h z 4 . For brevity, write z = Z(M\T; 7). 

The rank function of a matroid is submodular, so — rjn is supermodular 
and the probability distribution fi (defined by fi(A) = q~ r - M ^'y(A)/Z(M\T, 7)) 
associated with the random cluster model satisfies the condition for Fortuin- 
Kasteleyn-Ginibre (FKG) inequality [2J Section 19, Theorem 5], provided q > 1. 
Let f — 1 — e(A, {pi}) and g — 1 — e(A, {^2})- Note that these are both monotoni- 
cally increasing, as A grows. (This can either be seen directly, or as a consequence 
of submodularity of r^.) So the quantities / and g are positively correlated, i.e., 

Px.ifg = 1) = E^(fg) > E,(f) E,(g) = P T(t (f = 1) Pr^g = 1), 

which is equivalent to 

Zq . Zq + Z\ Zq + Z 2 

— > : x , 

z z z 

which in turn is equivalent to 

(23) z z > (z + Zi){zq + z 2 ). 



ISING MODEL ON A REGULAR MATROID 17 

Consider the first inequality we are required to establish, namely Zqz^ > z x z 2 . 
Using the linear transformation (I2"0"j) . we may express this as an equivalent inequal- 
ity in terms of 

z(Az + 2(zi + z 2 + z 3 ) + z 4 ) > (2z + 2zi + z 2 + z 3 + z 4 )(2z + z x + 2z 2 + z 3 + z 4 ), 

where the four linear factors may be read off from the appropriate rows (first, last, 
second and third, respectively) of the matrix V. Applying the definition of z we 
obtain the equivalent inequality 

z(z + 3z + Zi + z 2 + z 3 ) > (z + Zq + Zi)(z + z + z 2 ), 

which further simplifies, through cancellation, to 

(24) z z + z 3 z > (z + Zi)(zo + z 2 ). 

Now (J23|l implies fl24l . and we are done, since the other two advertised inequalities 
follow by symmetry. □ 

Denote by X 3 the 6-element matroid with ground set {pi,P2,P3, ei, e 2 , e 3 }, whose 
circuit space is generated by the circuits {pi, ei}, {p 2 , e 2 }, {p 3 , e 3 }, and {^1,^2,^3}- 
The matroid X 3 can be thought of as the cycle matroid of a certain graph, namely, 
the graph with parallel pairs of edges {pi, ex}, {p 2 , e 2 } and {p 3 , e 3 } in which edges 
Pi, p 2 and p 3 form a length-3 cycle in the graph. 

Let T = {pi,P2iVz\- We start by showing that, as long as a signature (s x , s 2 , s 3 ) 
satisfies certain equations, which Lemmas [8] and [9] will guarantee, then it is straight- 
forward to compute a weighting d so that the weighted matroid (X3, 8) has signa- 
ture cr(X 3 ;T,7) = (s 1 ,s 2 ,s 3 ). 

Lemma 10. Suppose Si, s 2 and s 3 satisfy 

(25) 2 + si - s 2 - s 3 > 0, 2 - si + s 2 - s 3 > 0, 2 - s 1 - s 2 + s 3 > 0, 

(26) si + s 2 + s 3 - 3 > 0, 
(27) 

Si+s 2 +s 3 -s 2 s 3 -2 > 0, S1 + S2 + S3-S1S3-2 > 0, and Si + s 2 + s 3 -s 1 s 2 -2 > 0; 

then there are non-negative weights d 1; d 2 and d 3 such that, for any weight function 
6 with 5 ei = di, 5 e2 = d 2 and S e3 = d 3 , a(Z 3 ;T,5) = (si,s 2 ,s 3 ). The values d\, d 2 
and d 3 can be computed from s\, s 2 and s 3 . 

Proof. Define 

51 = 2 + si - s 2 - s 3 

5 2 = 2 — Si + s 2 — S3 

5 3 = 2 — si — s 2 + s 3 
R = si + s 2 + s 3 - 2. 
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Define the weights di, d 2 , d 3 for 8 as follows: 

di = -1 + y/RS 1 /S 2 S 3 

d 2 = -l + V /RS 2 /S 1 S 3 

d 3 = -l + ^/RSz/SA, 

By inequalities (j2"5|) and (]2"fi]). i?, Si, and S 3 are all strictly positive. Thus 
di,d 2 ,d 3 are well defined. Finally 

RS 1 - S 2 S 3 = 4{ Sl + s 2 + s 3 - s 2 s 3 - 2) > 0, 

where the inequality is (|27p . Thus di, and hence, by symmetry, d 2 and d 3 , are all 
non-negative. 

Let Y = d\d 2 + d\d 3 + d 2 d 3 + d\d 2 d 3) and note that 

Z(X 3 \T; <5) = q-° + ^(di + d 2 + d 3 ) + q~% and 

Z(X 3 /pi\p 2 ,p 3 ; (5) = q-°(l + di) + q-\d 2 + d 3 + F). 

Substituting for d±,d 2 ,d 3 in these expressions, using the helpful identity 

Y + d 1 + d 2 + d 3 + 1 = (1 + di)(l + d 2 )(l + d 3 ) = Ry/R/S&Sa, 

we obtain 

Z(X 3 \T; 5) = + ^ + S 2 + S 3 ) ^R/S 1 S 2 S 3 = V /R/S 1 S 2 S 3 , 

and 

Z(l 3 / Pl \p 2 ,p 3 ; 8) = l(R + S 1 ) V /R/S 1 S 2 S 3 = s lV /R/S 1 S 2 S 3 , 

and hence Z(I 3 /p 1 \p 2 ,p 3 ; 5) / Z(Z 3 \T; d) = Si. By symmetry, similar identities 
hold for the other components of the signature of (X 3 , <5). Summarising, a(X 3 ; T, 5) = 
(si,s 2 ,s 3 ), as desired. □ 

Temporarily leaving aside the issue of approximation, the way that we will use 
Lemma [10] is captured in the following corollary, which is analogous to Corollary [7J 

Corollary 11. Suppose that (-Mi,7i) and (.M 2 ,7 2 ) are weighted binary matroids 
with E(M.i) n E(A4 2 ) = T = {pi,p 2 ,p 3 } , and suppose also that T is a circuit in 
both .Mi and M. 2 . Let (si,s 2 ,s 3 ) = o-(A4 2 ;T,j 2 ) ■ Then there are non-negative 
weights d\, d 2 and d 3 such that, for any weight function 5 with 5 ei = d\, 5 e2 = d 2 
and 5 e3 = d 3 , 

Z{M 1 A M 2 ; 7i A 7a ) = (Z{M X A X 3 ; j 1 A S), 

where 

C = Z(M 2 \T;j 2 )/Z(l 3 \T;S) = Z(M 2 \T;- f2 )/ y /R/S 1 S 2 S 3 , 

in the notation of the proof of LemmaUU The values d\, d 2 and d 3 can be computed 
from si, s 2 and s 3 — they do not otherwise depend upon (.Mi,7i) or (.M 2 ,7 2 ). 
Moreover, the values R, Si, S 2 and S 3 are byproducts of this computation. 
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Proof. As before, for i = 1,2, let 

z* = (Z(Mi\T; 7J, Z(Mi/pi\p2,P3] 7i ), Z(Mi/p 2 \pi,p 3 ; 7i), Z(M i /p 3 \p 1 ,p 2 ] 7i)) T 

Let s = (1, si, s 2 , s 3 ) T . By Lemma El 

Z(.Mi A M 2 ; 7i A 7a ) = (z^Dz 2 = Z(M 2 \T; >y 2 ){z l ) T Ds. 

Similarly, if ct(X 3 ; T, <5) = (si, S2, S3), then 

Z(Mi A X 3 ; 71 A <5) = Z(X 3 \T; <5)(/) T £> S . 

To simplify the notation (and avoid confusing the index "2" in z 2 with an expo- 
nent), let z denote the vector z 2 . Now, from inequalities (hi) and (iv) of Lemma El 
we see that 2{zq-\-z\) > z± + z 2 + z 3 , so Equation (|25|) is satisfied. By inequality (ii) 
of Lemma El Equation ( 126]) is satisfied. Finally, from Lemma [9] and identity (iv) 
of Lemma El 22^3 < zq{z\ + z 2 + z 3 — 2zq), so Equation ( 127]) is satisfied. So, by 
Lemma [TO] the weights d±, d 2 and d 3 can be computed and cr(Z 3 ; T, S) = (si, S2, S3). 
The value of Z(X 3 \T; S) is calculated in the proof of Lemma [TOl □ 

The problem with using Corollary [TTJ to replace the complicated expression 
Z(JAi A A4 2 ; 7 1 A 7 2 ) with the simpler Z(Aii A X 3 ; 7l A (5) is that, in general, 
we will not be able to compute the necessary values si, s 2 and s 3 . Instead, we will 
use our FPRAS recursively to approximate these values. Thus, we need a version 
of Corollary [TT] that accommodates some approximation error. Unfortunately, this 
creates some technical complexities. We will use the following lemma. 

Lemma 12. Suppose that (.Mi, 7l ) and (M. 2 ,~f 2 ) are weighted binary matroids 
with E(A4i) H E(M. 2 ) = T = {pi,p 2 ,p 3 } , and suppose also that T is a circuit in 
both Mi and M 2 . 

Suppose that e < 1 and that g is a sufficiently small positive constant (g = 
1/6000 will do). Suppose that, coordinatewise, 

(28) e' £S a(M 2 ; T, 7a ) < {s u §2, §3) < e £S a(M 2 ; T, 72 ). 

Then there are non-negative weights d\, d 2 and d 3 such that, for any weight func- 
tion d with 5 ei = d\, S e2 = d 2 and 5 e3 = d 3 , 

e~ e Z{M 1 A M 2 ; 7l A 7a ) < C^(A^i A X 3 ; 7l A5)< e £ Z{M 1 A M 2 ; 7i A 72), 

where 

C = Z(M 2 \T ]l2 )/ '^R/S^S 3} 
in the notation of the proof of LemmalTU. The values d\, d 2 and d 3 can be computed 
from §i, s 2 and s 3 — they do not otherwise depend upon (A^i, 7l ) or (Ai 2 ,j 2 )- 
Moreover, the values R, S\, S 2 and S 3 are byproducts of this computation. 

Proof. Let (7*1, r 2 , r 3 ) = a(Ai 2 ; T, 7 2 ). First, using Lemma [TBI in the appendix with 
X = sg, we can use §1, s 2 and s 3 to compute si, S2 and s 3 satisfying 

(29) e'^Si <r { < e 66£S s h 
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and Equations (I2"5l) . (j2B|) and (I2"T|) . To see that Lemma [TBI applies, note that Equa- 
tion (jinj) follows from Lemma E] (ii). Equation (jlT]) is analogous to Equation f [2"5"j) 
and is established in the same way as Equation ( I2"5"j) is established in the proof 
of Corollary [TTJ Similarly, Equation (138]) is analogous to Equation (|27|) and is 
established as in the proof of Corollary [TTJ (It is not clear that t\ < r 2 < r 3 , 
but this can be thought of as a renaming inside Lemma [16] Furthermore, the 
corresponding assumption si < §2 < S3 is without loss of generality, since, if the 
assumption does not hold, then these values can be swapped without violating the 
proximity of to ?y) 

By Lemma [TU| si, s 2 and s 3 can be used to compute non- negative weights di, d 2 
and d 3 such that, for any weight function d with S ei = d%, 5 e2 = d 2 and 5 e3 = d 3 , 
cr(X 3 ; T, 6) = s 2 , s 3 ). As in the proof of Corollary [TTJ this implies 

A X 3 ; 71 A «5) = Z(X 3 \T; ^(^flte, 

where s = (1, s 1; s 2 , s 3 ) T and Z(1 3 \T; S) = y/R/S&Ss for some byproducts R, 
Si, S 2 and S3 of the computation. 

Now let r = (1, r±, r 2 , r 3 ) T . Then, as in the proof of Corollary [TTJ 

Z(Mi A M 2 ; ii A 72 ) = Z(M 2 \T; 72 )(z 1 ) T Dr. 

Also, since 66g < 1/56, e _£//56 Sj < 7~j < e £//56 Sj. The result then follows from 
Lemma [T5] in the appendix, since z, s and r have positive entries, and satisfy 
1 < Zi/z < 2, 1 < Si / s < 2, and 1 < ri/r < 2 for i G {1,2,3} by Lemma 1 (i) 
and (ii). □ 

We also need a lemma, similar to Lemma [12], that is appropriate for 2-sum. 

Lemma 13. Suppose (A^i, 7l ) and (A4 2 , 72 ) are weighted binary matroids with 
E(A4i) Pi E(Ai 2 ) = {p}, where p is not a loop in either M. x or Ai 2 - Let z = 
Z(A / l 2 \p;7 2 ) and let z\ = Z(Ai 2 /p', 7 2 )- Suppose that e < 1 and that g is a 
sufficiently small positive constant ( as in Lemma UM — here it suffices to take 
Q < 1/160J. Suppose that e~ ee Zi < z% < e £e Zi fori e {0, 1}. Then there is a weight 
d > stzc/i £/&a£ ; /or even/ weight function 6 with 5 e = d, 

e- £ Z(M 1 A M 2 ; 7l A 7a ) < A X 2 ; 7l A 5) < e^A^ A A^ 2 ; 7l A 72 ), 

where ( = 2(2 + d) _1 zo- X/ie fa/we d can be computed from zq and z\. It does not 
otherwise depend upon (A^i, 7l ) or (A^ 2 , 72 ). 

Proof. This is similar to the proof of Lemma [T2] but easier. First, observe that 
zq > and 1 < zi/zq < 2 (Equation (TH|) ). The first step is to use z and Zi to 
compute Zq and z[ such that Zq > and 1 < z[/z' Q < 2 and e^ A£Q z' i < Zi < e i£Q z' i . 
This is straightforward. Just set z' = z e ee . If z\ < z' then set z[ = z' . If z\ > 2z' Q 
then set z[ = 2z' . Otherwise, set z[ = z x . (Note that g > e~ 2es ^ > e' 2ee , so if 
zi < Zq then zi = e^z < e ee e 2£e zi < e 4 ^. Similarly, | < e 2ee ^ < 2e 2ee , so if 
Z\ > 2zq then z[ = 2e ee z > e e£, e _2e£l Zi > e~ 2e£l z 1 .) As in the proof of Lemma [HI let 
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d = 2(z[ — z'q) I (2z' — z[). Let s' = z[/ z' and note that, for every weight function 6 
with 5 e = d, ~^ 2 ^ P ' S J = s'. Let s = z\/zq. As in the proof of Corollary [7J note that 

Zj \J-2 \p,0 ) 

Z{M, A M 2 ; 7l A 7a ) = (Z{M x \p; 7l ), Z{M x /p; 7l )) Q H z , 

and 

Z(^ 1 AX 2 ; 7l A<5) = (Z(Al 1 \p;7 1 ),Z(Al 1 /p;7 1 ))r_ 2 1 fV)(l + f)- 

Since e~ 8ee s' < s < e 8ee s', similar to the proof of Lemma [T5l 
e-^(Z(M 1 \pn 1 ),Z(M 1 /p;ji)) (_\ (j) 

□ 

4.1. Simple Delta-sums. We complete this section by investigating the (simple) 
way in which the special matroids I 2 and X3 interact with delta-sums with |T| = 1 
and \T\ = 3, respectively 

Suppose (M.j'y) is a weighted binary matroid on a ground set E(M.) of m el- 
ements, and with distinguished element p, which is not a loop. Consider the 
delta-sum Ai A X 2 . The ground set of this matroid also has m elements, and it 
shares all but one element with E(A4). Thus, we have a natural correspondence 
between the ground sets of the two matroids. We claim that under this correspon- 
dence the two matroids are the same, and for this it is enough to verify that they 
have the same circuit space. Note that {p, e} is the unique non-empty cycle in I 2 - 
For any cycle C in A4, exactly one of C = C or C = C © {p, e} is a cycle in 
A4 AX 2 . The mapping C \- > C is invertible, and is the required bijection between 
cycles in M. and those in M. AX 2 . Now suppose 8 e = d, and and let j' be derived 
from 7 by assigning 7 P = d. Then Z(Ai A X 2 , 7 A 5) = Z(J\A,~/'). 

A similar observation applies to X 3 under delta-sum with \T\ = 3. Suppose 
(A4,7) is a weighted binary matroid with distinguished elements T = {pi,P2,P3}, 
where T is a cycle in Ai. Consider the delta-sum Ai A X3. As before, there 
is a natural correspondence between the ground sets of the two matroids. Let 
Ci,C 2 , C3 be the three 2-circuits in X 3 including elements Pi,P2,P3, respectively. 
Any cycle C in A4 can be transformed in a unique way to a cycle C in Ai AX3, by 
adding a subset of circuits from {Ci, C 2 , C 3 }. The mapping C 1— > C is invertible, 
and is a bijection between cycles in Ai and those in Ai AX 3 . Now let 7' be derived 
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Step 1: If Ai is graphic, cographic or R w then estimate Z(Ai,^) directly. 

Step 2: Otherwise use Seymour's decomposition algorithm to express Ai as 
Mi A Ai 2 , where A is a 1-, 2- or 3-sum. Recall that E{M) = E{M\) © 
E{Ai 2 ). Let T = E{M 1 )nE{M 2 ), E l = E{Mx)-T and E 2 = E{Ai 2 )-T. 
Noting E(M) = Ex U £ 2 , let 7 X : #i -»■ M+ and 7 2 : £ 2 E + be the 
restrictions of 7 to i^i and £"2- Assume without loss of generality that 
\E(Ai 2 )\ < \E(Aii)\ (otherwise swap their names). 

Step 3: If A is s 3-sum then let T = {pi,p 2 ,Ps} (say). Execute Steps 4a- 
7a in Figured If A is a 2-sum then let T = {p}. Execute Steps 4b-7b 
in Figure [3l If A is a 1-sum then recursively estimate Z(A / fi;7i) with 
accuracy parameter e:|A1i|/|A1|, and Z(Ai 2 ; , y 2 ) with accuracy parameter 
^l-A^I/l-M I, and return the product of the two, which is an estimate of 
Z(M,~f). 



FIGURE 1. Algorithm for estimating the Ising partition function 
of a regular matroid Ai given accuracy parameter e < 1. g is a 
sufficiently small positive constant which does not depend upon Ai 
or e. See Lemmas fT2l and [131 

from 7 by assigning j pi = S ei , 7 P2 = 5 e2 and 7 P3 = 5 e3 . Then Z(A4 AX 3 , 7 A S) = 
Z(M.j'). 

5. The algorithm 

We now have all the ingredients for the algorithm for estimating Z(Ai;~f) = 
Z(A4;2,j), given a weighted regular matroid (Ai,j) and an accuracy parame- 
ter e. The base cases for this recursive algorithm are when Ai is graphic, cographic 
or R w . In these cases we estimate Z(Ai; 7) "directly", which means the following. 
If Ai is Rio then we evaluate Z(Ai;~f) by brute force. If Ai is graphic, we form 
the weighted graph (G, 7) whose (weighted) cycle matroid is (Ai,j). Then the 
partition function of the Ising model on (G, 7) may be estimated using the algo- 
rithm of Jerrum and Sinclair [8]. If Ai is cographic, then its dual Ai* is graphic, 
and 

Z{M; 1 ) = lEq~™ {E) Z{Ai*-n*)i 
where E = E(Ai), and 7* is the dual weighting given by 7* = q/^ e = 2/ , y e for 
all e G E(Ai) [131 4.14a]. (Ground set elements e with 7 e = do not cause any 
problems, because they can just be deleted.) Then we proceed as before, but using 
(Ai*,"y*) in place of (Aij'y). The proposed algorithm is presented as Figured] 

Using the guarantees from Lemmas fT2l and [TBI it is easy to see that the algorithm 
is correct. That is, given a regular matroid Ai and an accuracy parameter e < 1, 
the algorithm returns an estimate Z satisfying e~ e Z(Ai,~y) < Z < e e Z(Ai,~y). 
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Step 4a: Recusively estimate zq = Z(Ai 2 \T;~f 2 ), z i = Z(Ai 2 /pi\p 2 ,p3] 72X 

z 2 = Z(M 2 /p 2 \p 1 ,p 3 ;^ 2 ) and z s = ^(MzMVi,^; 72) witn accuracy 
parameter e q\M2\ / (^\M\) . 
Step 5a: Using Lemma fT2| compute d±, d 2 and such that, for any weight 
function 8 with 5 ei = d%, 5 e2 = d 2 and 5 e3 = d 3 , 

e -c\M 2 \/(2\M\)z(M;~r) < CZ(M 1 Al 3 ;~ fl AS) < e £ \ M ^^ M ^Z(M^). 

Note that our estimate for z gives an estimate for ( = z /\/R/ S\S 2 Ss 
with accuracy parameter at most e|A^2|/(2|A^|). (R, Si, S 2 and 5*3 are 
byproducts of the computation of d\, d 2 and d 3 .) 
Step 6a: Recall from Section O that Z{M X AZ^7j A 8) = Z(M u j'), 
where 7' is derived from j 1 by assigning j' pi = 6 ei , 7p 2 = 5 e2 and 7p 3 = 5 e3 . 

Step 7a: Recursively estimate Z(A4i;7') with accuracy parameter e(|.M| — 
|A^2|)/|A^| and multiply it by the estimate for ( from Step 5a. Return this 
value, which is an estimate of Z(A4,j). 



Figure 2. The 3-sum case. 



Step 4b: Recursively estimate z = Z(.M2\p; 72) an d z\ = Z(A4 2 /p;j 2 ) 

with accuracy parameter (e^|A^2|)/(2|A / l|). 
Step 5b: Using Lemma [T3| compute d such that for every weight function 8 

with 5 e = d, 

e -e\M a \/{2\M\)z( M .^ < (z{M x AZ 2 ; 7l AS) < e elM2l/(2lMl) Z(M; 7). 

Note that our estimate for zq gives an estimate for £ = 2(2 + d)~ 1 zo with 
accuracy parameter at most (e|A^2|)/(2|A^|). 
Step 6b: Recall from Section O that Z(Mi AI 2 , 7l AS) — Z(M%,j'), 
where 7' is derived from "y 1 by assigning 7^ = d. 

Step 7b: Recursively estimate Z{M.\\~i') with accuracy parameter e(|.M| — 
|A^2|)/|-^t| an d multiply it by the estimate for £ from Step 5b. Return this 
value, which is an estimate of Z(Ai,~y). 



Figure 3. The 2-sum case. 

The rest of this section shows that the running time is at most a polynomial in 
\E(M)\ and £~\ 

Let Tdecomp(^) = 0(m Odccomp ) be the time complexity of performing the Seymour 
decomposition of an m-element matroid (this is our initial splitting step), and let 
Tb ase (m, e) = 0(m abasc e~ 2 ) be the time complexity of estimating the Ising partition 
function of an m-edge graph. (From [HI Theorem 5], and the remark following it, 
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we may take abase = 15.) Denote by T(m, e) the time-complexity of the algorithm 
of Figure [TJ Recall that g is a sufficiently small positive constant which does not 
depend upon M. or e and is described in Lemmas d2] and [T3J we can take it to 
be q = 1/6000. The recurrence governing T(m,e) is now presented, immediately 
followed by an explanation of its various components. 



T(m, e) < T decomp (m) + max |T base (m, e) 



max (T(m - k + 3, + 4T(k, . 

4<k<m/2 



e(m-fc-l) ^ | 2T( U ggfetll ' 

(T(m-fc,^)+T(fc,S))" 



max (T(m - fc + 1, **=*=2) + 2T(*, ^)) , 

2<fe<m/2 



max 

l<fc<m/2 

The four expressions within the outer maximisation correspond to the direct case, 
the 3-sum case, the 2-sum case, and the 1-sum case, respectively. The variable k is 
to be interpreted as the number of ground set elements in Ai that come from A4 2 
(and hence m — k is the number that come from .Mi). Thus, in the case of a 3-sum, 
for example, \E(Mi)\ = m - k + 3 and \E(M 2 ) \ = k + 3. Note that Step 2 of the 
algorithm ensures k < m/2. The lower bounds on k come from the corresponding 
lower bounds on the size of matroids occurring in 3-sums, 2-sums and 1-sums. 

We will demonstrate that T(m, e) = 0(m a e~ 2 ), where a = max{«b a se, adccomp + 
1,43}. Specifically, we will show, by induction on m, that T(m,e) < Cm a e~ 2 , for 
some constant C and all sufficiently large m. In the analysis that follows we do 
not attempt to obtain the best possible exponent a for the running time, instead 
preferring to simplify the analysis as much as possible. It would certainly be 
possible to reduce the constant 43 appearing in the formula for the exponent, but 
there seems little point in doing so, as the best existing value for abase is already 
too large to make the algorithm feasible in practice. 

Note that by choosing C sufficiently large in the time-bound T(m, e) < Cm a e~ 2 , 
we may ensure that this bound on time-complexity bound holds for any m from 
a finite initial segment of the positive integers, specifically for m G {1,2,..., 19}. 
For the inductive step, substitute T(m,e) < Cm a e~ 2 into the right-hand side of 
the above recurrence. We need to verify 

(30) Cm a e~ 2 > T dccomp (m) + max |T base (m, e), 

(31) 4 max^ (C(m - k + ^{j^r^) 2 + 4^(^) 2 ), 

(32) a max /2 (C(m - k + l) Q (^Ty) 2 + 2Ck"(^) 2 ) 

(33) ^max 2 (C(m - kr(^) 2 + Ck«{£) 2 ) }, 
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for m > 20, as this will imply T(m,e) < Cm a e~ 2 , completing the induction step. 
Effectively, there are four independent inequalities to verify, numbered f[3"U|) - f[3"3"j) . 
The first inequality is immediate (for a sufficiently large constant C), since we are 
assuming a > max{a dccomp , a ba sc}- 

The second inequality, namely (13T|) . requires some work, but the remaining two 
will then follow easily. First, note the following simple estimate: 

m-k + 3 m/2 + 3 12 20 

34 -—— < ' = 1 + < 1 + — 

m — k — 3 m/2 — 3 m — 6 m 

where we have assumed that k < m/2 and m > 20. The inequality we wish to 
establish is equivalent, under rearrangement, to 

Cm«e- 2 -C(m-k+ ?>T{ 7{ ^) 2 - ^{j^f > T dccomp (m), 

for all k with 4 < k < m/2. Noting fl34l) . it is enough to show 

(35) Cm 2 e~ 2 [m a ~ 2 - (m - k + 3) a ~ 2 (l + 20/m) 2 - QAk a ~ 2 g- 2 ] > T dccomp (m), 

for all k with 4 < k < m/2. Regarding the left-hand side of fl35|) as a continuous 
function of a real variable k, and taking the second derivative with respect to k, 
we that the left-hand side is a concave function of k. It is enough, then, to check 
that (J31>D holds at k = 4 and k = m/2. 

When k = 4, the inequality follows from the following sequence of inequalities: 

(36) Cm 2 e- 2 [m a ~ 2 - (m - 1) Q ~ 2 (1 + 20/m) 2 - 64 x 4 a ~ V 2 ] 

> Cm 2 e- 2 [m a - 2 - (m - 1) Q - 2 (1 + 1/m) 40 - 4 a+ V 2 ] 

= Cm 2 e~ 2 [m a ~ 2 - m fl " 2 (l - 1 /m) a " 2 (l + 1 /m) 40 - 4 Q+1 ^ 2 ] 

(37) > Cm 2 e- 2 [m a ~ 2 - m a " 2 (l - 1/m) - 4 a+ V 2 ] 
= Cm 2 s- 2 [m a - 3 -A a+l g- 2 ] 

(38) > ±C7m a - V 2 

(39) > T dccomp (m), 

where inequality fl37j) is a consequence of a > 43, inequality (1391 of a > a dccomp -|-l, 
and inequality fl38|) comes from a comparison of the two terms, noting a > 43, 
q > 1/6000 and m > 20. 

When k = m/2, inequality fl35|) is established as follows: 

(40) Cm 2 e' 2 [m a ~ 2 - (m/2 + 3)°- 2 (l + 20/m) 2 - (m/2) a - 2 (8/ g) 2 ] 

= Cm a e- 2 [1 - (i + 3/m) a - 2 (l + 20/m) 2 - (l/2) a - 2 (8/ g) 2 ] 

(41) > Cm a e~ 2 [1 - 10" 7 - 10~ 2 ] 
> \Cm a e- 2 

(42) > T dccomp (m), 
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where (Jl"2"]l is a consequence of a > «decomp, and flH} of a > 43, £ > 1/6000 and 
m > 20. 

The above calculations may easily be adapted to cover inequalities (}3"2"]) and ( 133j) 
for the cases of 2-sum and 1-sum. Since estimate (1311 applies as well to (m — + 
l)/(m — k — 1), the analogue of (|35|) in the 2-sum case is 

Cm 2 e" 2 [m"- 2 - (m - k + l) a ~ 2 (l + 20/m) 2 - 8k a - 2 g- 2 ] > T dccomp (m). 

Again, we need only verify this at the extreme values of k, namely k = 2 and 
k = m/2. The specialisation to k = 2, 

CmV 2 [m a - 2 - (m - 1) Q ~ 2 (1 + 20/m) 2 - 2 Q+ V 2 ] > T dccomp (m), 

can be seen by comparison with (135"j) . and the one for k = m/2, 

Cm 2 e- 2 [m a ~ 2 - (m/2 + l) a ~ 2 (l + 20/m) 2 - 8(m/2) Q ~V 2 ] > T dccomp (m). 

by comparison with ( 140]) . 

Finally, the analogue of (1351) in the 1-sum case is 

Cm 2 e- 2 [m a - 2 - (m - k) a ~~ 2 - k a ~ 2 ] > T dccomp (m), 

which again needs to be verified at the extreme values k = 1 and k = m/2. As 
before, these can be seen by comparison with (136]) and (T4T1 . 
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6. Appendix 

6.1. Part one. The first part of the appendix contains some additional calcula- 
tions needed in the proofs of Lemmas H] and [5] to provide lower bounds on the 
quantity c(A, B). We believe that several people have already done this calcula- 
tion. Essentially, it shows that the two special cases of delta-sum that we consider 
can be constructed via a "generalized parallel connection" , which is a special amal- 
gam of M.\ and M.2 called a "proper" amalgam. This fact is well-known. See, for 
example [TP"} Section 12.4] or pQ. We couldn't find the calculation written down 
anywhere, and our result relies on it, so we include the calculation here. We use 
some ideas from [T2J Section 2]. 

6.1.1. Lemma gj- iV = A/i- Recall that c(A, B) = r M (A U B) - r Ml (A) - r M2 (B). 
Our goal is to show 

> f-1 iff ei {A, M) = and e 2 (B, {p}) = 0; 
— [ otherwise. 

Our basic strategy will be to identify sets A" C A and B" C B such that 
r Ml (A") = \A"\ and r M2 (B") = \B"\ and 

(43) e 1 (A",{p}) + e 2 (B",{p})>0 

Equation flU implies that r M (A"UB") = \A"UB"\. (Otherwise, A"UB" contains 
a non-trivial cycle C. By the definition of delta-sum, it is of the form C = C\ © C 2 
with Ci G C(Mi) and C 2 G C(M 2 ). Since A" is independent in Mi and B" is 
independent in M. 2l neither C\ nor C 2 is trivial. Now C\ is a subset of A" U {p}, 
so A" U {p} is dependent in .Mi, but A" is independent in A^i, so ei(A",p) = 0. 
Similarly, e 2 (B",p) = 0, which is a contradiction.) Then 

c(A, B) =r M (AUB)- r Ml (A) - r Mi (B) 

(44) > \A"\ + \B"\ - r Ml (A) ~ r M2 {B). 

To implement the strategy, we first choose B" C B with r^ 2 (i?) = r^ 2 (S") = 
\B"\. Choose A' C A with r Ml (A) = r Ml (A') = \A'\. If e 1 (A', {p}) + e 2 (B", {p}) > 
(which is certainly true as long as ei(A, {p}) + e 2 (B, {p}) > ) then we take 
A" = A' so (HU) gives c(A, B)>0. 
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Otherwise, A' U {p} contains a circuit. Let e / j) be a member of this circuit 
and take A" = A' - e. Then r Ml (A") = \A"\ = \A'\ - 1 and r Ml (A" U {p}) = \A'\ 
so ei(A", {p}) > 0. Equation (JSj) gives Then c(A, 5) > -1. 

6.1.2. LemmaS AT = A3. Our goal is to show 

( q 2 q q q l\ 

q q 1 1 1 

C < q 1 q 1 1 

g 1 1 q 1 

Vi 1 1 1 v 

where (7^ = q~ c ( A ' B \ for < j, k < 4, where ACEj and B C E 2 are any sets 
satisfying <£>j(A) and ipl(B). 

The basic strategy is similar to the one used in Appendix 16.1.11 We wish to 
identify sets A" C A and B" C 5 such that r^A") = and r M2 (B") = \B"\ 
and, for every p,,- G T, 

(45) e 1 (A",{p i }) + e 2 (5",{p i })>0. 

Equation (T4"5]) implies that ^(^'Ui?") = |A"US"|. The reasoning is a little more 
complicated than in Appendix 16. but it is similar. If A"UB" is not independent 
in M. then, as before, A" U B" contains a non-trivial cycle C = C\ © C2 where 
the non-trivial cycle C\ G C(Mi) is a subset of A" U T and the non-trivial cycle 
C 2 G C(A^ 2 ) is a subset of 5" U T. Note that C = (C l © T) © (C 2 © T) so we 
can assume, without loss of generality, that \C\ D T| < 1 (otherwise, replace C\ 
with Ci © T and replace C 2 with C 2 © T). Also, Ci fl T is nonempty. Otherwise, 
C\ C A", contradicting the independence of A" in .Mi. So suppose C\ flT = {pj}. 
We can deduce that C7 2 fl T = {p^}, since C avoids T. Now, as in Appendix 16. 1.11 
we get a contradiction. Since C\ C A" U {p^}, but A" is independent in A^i, 
we get ei(yl",pj) = 0. Similarly, e 2 (5",p J ) = 0, which is a contradiction. As in 
Appendix 16.1.11 implementing the basic strategy allows us to conclude 

(46) c(A,B) > \A"\ + \B"\-r Ml (A)-r M2 (B). 

To implement the strategy, we choose B' C B with r M2 (B) = r M2 (B') = \B'\ 
and we choose A' C A with rj^ x (A) = r^ 1 ( J A / ) = \A'\. 

If, for all pj G T, ex(A', {pj}) + e 2 (B',{pj}) > 0, (which is certainly true as 
long as, for all pj G T, ei(A, {pj}) + e 2 (S, {pj}) > ) then we take A" = A' and 
B" = B' so fT46|) gives c(A, -B) > 0. This accounts for all of the "1" entries in the 
matrix above. 

Suppose, instead, that e\(A', {pi})+e 2 (B' ', {pi}) = but e\(A\ {p 2 })+e 2 (B' ', {p 2 }) > 
and e\{A', {p 3 })+e 2 (i? / , {P3}) > 0. (This case arises for all cases with ei(A, {p 2 }) + 
e 2 (-B, {p 2 }) > and ei(A, {p 3 }) + e 2 (B,{p 3 }) > that we have not already 
covered.) In this case, A' U {pi} contains a circuit. Let e 7^ pi be a member 
of this circuit and take A" = A' - e. Then r Ml (A") = \A"\ = \A'\ - 1 and 
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r Ml (A" U {px}) = \A'\ so ex{A", {pi}) > 0. Taking 5" = B', Equation (JH gives 
Then c(A, B) > — 1. By symmetry (swapping pi with either p 2 or p 3 ), this accounts 
for all of the "g" entries in the matrix above. 

For the final case, suppose that, for all pj £ T, ei(A', {pj}) + ^(B', {pj}) = 0. 
Define A" as above so that ei(A", {pi}) > 0. As we note just after enumerating 
the possibilities P0, . . . , P4 in the proof of Lemma [5j it is not possible that exactly 
one of ei(A", {pi}), ei(A", {p 2 }), and ei(A", {^3}) is 1, so assume without loss of 
generality that ei(A", {p 2 }) > 0. Now B' U {p 3 } contains a circuit. Let / 7^ p 3 be 
a member of this circuit and take B" = B' — f so rj^ 2 {B") = \B"\ = \B'\ — 1 and 
e 2 (B", {P3}) > 0. Equation (H5T) now gives c(A, B) > —2, which accounts for the 
q 2 entry in the matrix. 

6.2. Part two. The second part of the appendix contains some technical lemmas 
needed in the proof of Lemma fT2l The lemmas are about approximation, and they 
don't add any intuition to the paper. 

We start with a well-known fact that we will use in the proof of both lemmas. 
(This follows directly from the series expansion of e.) 

Observation 14. If < e < 1 then l + £<e e <l + 2e. 

Our first lemma refers to the matrix D defined just before Lemma [5j 

Lemma 15. Suppose that z, s and r are column vectors in R 4 with positive entries 
satisfying 1 < z^/zq < 2, 1 < Si/so < 2, and 1 < r^/Vo < 2 for i £ {1,2,3}. 
Suppose that e~ £ Si < Vi < e e Si for some < e < 1. Then e~ 56e z T Ds < z T Dr < 
e 56e z T Ds. 

Proof. For any positive column vector v in R 4 satisying 1 < V{/vq, 

z T Dv = z v + (zi - z ){vi - v ) + (z 2 - z )(v 2 - v ) + (z 3 - z )(v 3 - v ) > Z V . 

Also, summing the absolute values of all of the monomials in z T Dv, and using 
Zi/z < 2 and Vi/v < 2, we get \ D i,j\ z i v j < 28z v . Then 

z T Ds — z T Dr = DijZi(sj — rf) 
hi 

< 2e^2 \ D i,j\ z i r i 

< 56e2;oro 

< 56ez T Dr, 

so z T Ds < (1 + 56e)z T Dr < e 56e z T Dr. □ 
Our second lemma refers to the equations in Lemma (TDl 
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Lemma 16. Suppose that r\ < r 2 < r 3 satisfy the following equations. 

(47) 2 + n - r 2 - r 3 > 0, 

(48) n + r 2 + r 3 - r 2 r 3 - 2 > 0, 

(49) 1 < n < 2, fori e {1,2,3}. 

Suppose that we are given values §\ < s 2 < s 3 satisfying e~ x ri < s~i < e x ri for 
i £ {1,2, 3} ; where x is a sufficiently small positive constant. Using s~i, § 2 and s 3 , 
we can compute values si, s 2 and s 3 satisfying 1 < Sj < 2 ; 

e" 66 ^ < r-i < e 66x Si, 

and Equations / f!3]) . ( f^) and 
Proof. Let 5 = 4e%. Note that 

(50) | S ._ r .|< ( j. 

For example, using Observation dU — Sj < 2x5j < 2xe x rj. Since < 2 (by 
Equation ( T49]) ) and x < 1, r,i — < 5. Similarly, — < 5. 
In each of two cases, we will compute S\, s 2 , and s 3 so that 

(51) 1 < si < s 2 < s 3 < 2, 

(52) 2 + Sl -s 2 -s 3 >0, 

(53) si + s 2 + s 3 - s 2 s 3 — 2 > 0, 

(54) \si-n\<66. 

Equations (1ST]) . (152]) and (153]) imply Equations 025]), (J2SJ) and (J27J). Also, Equa- 
tion (154]) implies e~ 66x Sj < < e 66x S; since Sj < rj+65 < rj(l+65) < r^e 6 " 5 < rje 66x 
and similarly r« < Sje 66x . 
Case 1: s 2 — si < 55. 

Take Si = s 2 = min(max(l, 5 2 ), 2 — 5) and s 3 = min(max(l, S 3 ), 2 — 5). ( 15T]) 
follows easily from the definitions since 2 — 5 > 1 and 5 2 < § 3 . (152]) and (153]) follow 
from the facts that s± = s 2 , Si > 1 and s 3 < 2 since, Si + si + s 3 — Sis 3 — 2 = 
(si - 1)(2 - s 3 ). To establish (jSJD note that 

sx < max(l, 5 2 ) < max(l, §i + 55) < max(l, r\ + 65) < r\ + 65. 

Also, since S 2 < r 2 + 5 < 2 + 5, 

si > s 2 - 25 > si - 25 > n - 35. 

Similarly, r 2 — 35 < s 2 < r 2 + 5 and r 3 — 35 < s 3 < r 3 + 5. 
Case 2: s 2 — s 1 > 55. 
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Take s± — §i + 45, s 2 = min(§ 2 , 2) and s 3 = mm(s 3 ,2). Equation f l5"Tj) follows 
since §i + 45 < s 2 — 5 < s 2 . flo^]) follows since 

2 + si - s 2 - s 3 > 2 + §i + 45 - s 2 - s 3 

>2 + ( ri -5)+45- (r 2 + 6) - (r 3 + 6) 

> 5 + 2 + r l - r 2 - r 3 

> 5. 

We now show that f)53p holds. Note that 

si + s 2 + s 3 - s 2 s 3 - 2 = si - (s 2 - l)(s 3 - 1) - 1, 

and the this quantity is increasing as a function of Si and decreasing as a function 
of s 2 and as a function of s 3 . Also, 

si = si + 45 > ri + 3(5, 

s 2 = min(s 2 , 2) < s 2 < r 2 + <5, and 

s 3 < ^3 + S. 

So 

si + s 2 + s 3 - s 2 s 3 - 2 > (n + 35) - (r 2 + 6 - l)(r 3 + 5 - 1) - 1 

= 35 - 5(r 2 - 1) - 5(r 3 - 1) - 5 2 + (n - (r 2 - l)(r 3 - 1) - 

> 35 - 5(r 2 - 1) - 5(r 3 - 1) - 5 2 , 

and this is at least 0, since r 2 — 1 < 1, r 3 — 1 < 1, and 5 < 1. 

To establish (154"]) . note that S\ = §i + 45 < r x + 55 and r x < si. Also, as noted 
above, s 2 — 5 < s 2 < s 2 and similarly s 3 — 5 < s 3 < s 3 . 

□ 
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